Generalisation of Gunion-Bertsch Formula for Soft Gluon Emission 
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We generalise the most extensively used Gunion-Bertsch formula for the soft gluon emission in 
a perturbative QCD. We show that the corrections arising due to this generalisation could be very 
important for the phenomenology of the hot and dense matter produced in the heavy-ion collisions. 
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The prime intention for ultra relativistic heavy-ion col- 
lisions is to study the behaviour of nuclear or hadronic 
matter at extreme conditions like very high temperatures 
and energy densities. A particular goal lies in the identifi- 
cation of a new state of matter formed in such collisions, 
the quark-gluon plasma (QGP), where the quarks and 
gluons are deliberated from the nucleons and move freely 
over an extended space-time region. Various measure- 
ments taken in CERN-SPS 1] and BNL-RHIC 043] do 
lead to 'wealth of information' for the formation of QGP 
through the hadronic final states. In the upcoming ex- 
periments at the CERN LHC, one is hoping to produce 
QGP during the first several fm/c of the collisions and 
to substantiate those evidences already found in the past 
as well in the recent experiments. 

Some quantitative key features of the plasma produced 
in such collisions include equilibration and its time, initial 
temperature, energy-loss and jet quenching of high ener- 
getic partons, and elliptic flow of hadrons and its scaling 
with the number of valence quarks. The Gunion-Bertsch 
(GB) formula @ for soft gluon emission has widely been 
used for various aspect of the heavy-ion phenomenology. 
To set the perspective we note that there are many pa- 
pers in the literature based on the GB formula. We 
recall: the sequence of events in hot glue scenario 0- 



12j, thermal equilibration 13 . 14$ . gluon chemical equi- 
libration [l3|, 0, Eij], parton matter viscosity [l?], EH, 
radiative energy-loss of high energy partons propagat- 
ing through a thermalised QGP [19rl21| etc., where the 
GB formula has extensively been used. The original GB 
formula was derived in Ref.0] for gluon emission from 
quark-quark scattering and later it was explicitly used in 
Ref.flH.l23 to derive the soft gluon emission from gluon- 
gluon scattering. Recently in Ref. 23j a correction to 
the GB formula for soft gluon emission was obtained. In 
present article we make an effort to generalise the GB for- 
mula for soft gluon emission from gluon-gluon scattering, 
i.e., gg — > ggg and find a more important correction than 
it is found in Ref. [23[ . We also show that in an appropri- 
ate limit the generalisation reduces to the GB formula. 
We further note that the results for similar inelastic pro- 
cesses can be obtained in a straightforward way by using 



our generalisation. 

For the process, gg — > ggg, there are 25 different 
topologies. We note that k\ and k 2 are momenta of the 
gluons in the entrance channel, k 3 and fc 4 are those for 
exit channel gluons whereas fcs is that of the emitted 
gluon. The invariant amplitude summed over all the fi- 
nal states and averaged over initial states for such process 
can elegantly be written [24[ as: 
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where N c is the number of color, g — y'Airas is the strong 
coupling, 

Af = (fci • k 2 ) 4 + (fci • k 3 ) 4 + (fa ■ fa) 4 + (fa • fa) 4 
+(k 2 ■ fa) 4 + (fa ■ fa) 4 + (k 2 ■ fa) 4 + (k 3 • fa) 4 
+ (k 3 ■ fa) 4 + (k 4 ■ fa) 4 , (2) 

V = (hi- k 2 )(h ■ k 3 )(h ■ k 4 )(h ■ k 5 ){k 2 ■ k 3 ) 

(k 2 ■ ki){k 2 ■ k 5 )(k 3 ■ k 4 )(k 3 ■ fc 5 )(fc 4 • k 5 ) , (3) 

and 

(ijklm) = (ki ■ kj)(kj ■ k k )(k k ■ k t )(ki ■ k m )(k m ■ ki) . (4) 

We now define the Mandelstram variables as 

s = {k x + k 2 f , t = (h - k 3 ) 2 , u={kx- fc 4 ) 2 , 
s' = {k 3 + k A f, t' = {k 2 -k i f, u' = {k 2 -k 3 ) 2 , 

and 

s s' t 

h-k 2 = -, fc 3 -fc 4 = — , h-k 3 = --, 

k 2 -ki = --, ki-ki = ~, k 2 -k 3 = -^-, 

fa.* = £±1+4, k 2 . k5 = ^^l, 

k 3 -k 5 = il±^, fa- fc5 = ^±±±^,(5) 
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with, 



s + t + u + s' + t' + u = 



Eq. ([T]) actually contains only twelve terms which are 
reduced from total 120 terms due to symmetry (24[. Us- 
ing ©, ^ and dU) the first two terms of (P) can be 
expressed as: 



2 y N 2 - 1 (ki ■ fc 3 )(*i • k 4 )(k 2 ■ k 4 )(k 2 ■ k 5 )(k 3 ■ h) 
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and, 



2 a N 2 -l(fc x - fcs)(*i • fe)(fc 2 • *4)(*a ' • fc 4 ) 
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After simplifying all the terms in this way, ([T]) can be 
written as 
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(8) 



t(s + u' + t')(s +u' + t) \ss' uu' 
and AT can also be obtained as 

N = ^(s i + t i + u i + s' 4 + t' 4 + u' 4 ) 

+ T6 \-( s + t + u ) 4 + ( s + *' + u ') 4 + ( s + *' + u ) 4 



+ (s + t + u') 4 



(9) 



Now for a soft gluon emission (k$ — > 0): i' — > t, s' — > 
s, it' — > it and we can express the transverse component 
of the momentum of the emitted gluon in the centre of 
momentum frame of k\ and k 2 as 



k\ = 



4(fei • fc 5 )(fc 2 • fc 5 ) (s + t + u)(s + u' + f) 



(s + t + u) 2 



(10) 



Using (IT0|) in ([8j one can obtain a complete expression 
for the three-body matrix element in terms of the two- 
body matrix element and a infrared factor as 
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where the full two-body matrix element ll|, l25j is given 
as 



9 A s 2 
2 9 ¥ 



and various coefficients are 
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Furthermore, eliminating u using (| 10[) and keeping 
terms upto C(l/fcj_) and 0(i 3 /s 3 ) one can write (TTT|) 
as 
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The above equation (IT4|) is a convergent series of t/s as 
can be seen below. It is now clearly decomposed in two 
factors: one is associated with the 2 — > 2 process used 
by Gunion and Bertsch 0] whereas the other one is the 
generalisation of the infrared factor for the emission of 
soft quanta. As we will see below, the first term in (fT4")) 

will lead to the GB formula 0] for the gluon multiplicity 

2 

distribution in an appropriate limit -5- w (FjT^Tp ' 
where t = qj_ >> k\ (q± is the transverse component of 
the momentum transfer). If the emitted gluo n is much 
softer than others it can then be regulated 15, 22| by the 
Debye screening mass, mr>- On the other hand the sec- 
ond, the third and the fourth terms, respectively, would 
correspond to the correction terms over the GB term in 
the same spirit. We also note that the first term can also 
be obtained by just using the scalar QCD approximation 
in 2 — > 3 process. 

If one considers only a\, 04, ciq, a~i in (jlip and substitute 
u = —a, one will end up with the expression derived in 
Ref. SI as 



12 3 2 |X„/ GB ^ (\ + -^), (16) 
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FIG. 1: The ratio R = 



7 



— as a function of \fs 

GB 



dt]dk ^ ' dr]dk 

at the center of momentum frame of gluon-gluon scattering 
at T = 200 MeV and a s = 0.3. 



FIG. 2: The ratio R = 
strong coupling, a s . 



drjdk-. 



dr)dk 2 . 



as a function of 



which misses the leading order correction term of 0(t/s) 
and also the right coefficient in the next order, 0(t 2 / s 2 ). 

Using (Q3]) and following Refs. [H, [22[, it is now very 
straight forward to obtain the soft gluon multiplicity dis- 
tribution in the midrapidity region as 
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where the GB formula can be obtained using ([15)) as 



dn„ 



dr/dk 2 ; 



Cao-s 



gb 



k 2 ± [(k ± - q A 



(18) 



"Di 



with the Casimir factor C 'a = 3 and the Debye screening 
mass mo = gT. We note that the first term in {JTTJ) 
would correspond to GB formula as obtained in Ref. [8j . 
It is also worth mentioning that this first term can be 
obtained by just using the scalar QCD approximation. 

Now the average value of the momentum transfer 
squared can be obtained as 
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where the maximum value of a 
along with s = (s) = 18T 2 . 

In Fig. [1] the ratio R 



is restricted by s > 4m 



drjdk' 



-I 



da ( . 



drjdk^ 



is plotted 



with the center of momentum energy, s/s of the gluon- 
gluon scattering for T = 200 MeV and a s = 0.3. With 



these values the lower limit of y/s ~ 0.8 GeV is restricted 
by the relation s > Am 2 D . As can be seen from Fig. [1] that 
(H3]) is a convergent series in t/s and it also indicates a 
significant improvement compared to that of Ref. [23j |. 
With the increase of s/s the correction terms decreases 
and the ratio approaches the GB formula for very large 
value of y/s. On the other hand the first and second order 
corrections are significant when \fs < 6 GeV. 

In Fig. ® the ratio R = Jfe-/ "'" 



drjdk'^ ' drjdk 2 ± 

a function of the strong coupling a 
maximum value of a 



is plotted as 

GB 

We note that the 



is restricted by s > Am 2 D . As can 



be seen the correction terms become important with the 
increase of a s . At a s — 0.3, the correction up to 0(t 3 /s 3 ) 
is around 50% over the usual GB formula. We have also 
compared our results with that of Ref. [23j , which shows 
a significant improvement in the range of a s displayed in 
Fig.H 

Fig. [3] displays the ratio R as a function of T in the 
units of the critical temperature T c . Here we have used 
the temperature dependent a s from Ref. [26| with two 
momentum scale Q = 2ttT (red curve) and AnT (green 
curve). We note that for the momentum scale Q = 2irT 
the value of the coupling, a s > 0.4 for T/Tc < 1.5, which 
is really a very high for any practical purposes. On the 
other hand for Q = 4irT, the values of coupling lie in 
the domain 0.2 < a s < 0.35 for the scaled temperature 
range, 1 < T/Tq < 6. Nonetheless, the value of a s and 
thus the lower bound oiT/Tc would again be restricted 
by the relation s > Am 2 D . As can be seen in Fig.[3l there 
is a sizeable contribution (> 40%) up to 0(t 3 /s 3 ) over 
the GB formula and a correction (> 30%) compared to 
Ref. (Hi in the temperature range 1 < T/T c < 5. This 
correction over the GB formula would be very important 
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FIG. 3: The ratio R = 



m. 



-i 



up to 0(t 3 /s 



displayed as a function of T/T c for temperature dependent ct 3 
with two momentum scales 2nT (red curve) and AnT (green 
curve). The corresponding results from Ref. are repre- 
sented by dotted (red) and dashed (green) curves. 



at the temperature domain relevant for the heavy-ion col- 
lisions. 

In summary, we have obtained a generalisation of the 
GB formula for soft gluon emission in a process gg — > ggg. 
We found that the correction terms are important at var- 
ious physical domains of temperature, coupling constant 
and the energy of gluon-gluon scattering. This generali- 
sation will be very apt for the phenomenology (viz., hot 
glue scenario, chemical equilibration of gluons, partonic 
matter viscosity, radiative energy-loss of energetic par- 
tons and jet quenching) of heavy-ion collisions and would 
improve the present understanding on various phenom- 
ena in this area. 
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